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Abstract 

A new solvable many-body problem of goldfish type is introduced and the 
behavior of its solutions is tersely discussed. 
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1 Introduction 

Notation 1.1. Hereafter N is (unless otherwise indicated) an arbitrary integer, 
N > 2, the (generally complex) numbers Zn = Zn (t) are the dependent vari¬ 
ables, t (“time”) is the independent variable, superimposed dots denote time- 
differentiations, and indices such as n, m, I run over the integers from 1 to A^ 
unless otherwise indicated (see for instance below in (la) the limitation non 
the N values of 1). Below we often omit to indicate explicitly the t-dependence 
of various quantities, when this can be done without causing misunderstandings. 
Hereafter N x N matrices are denoted by upper-case boldface letters (so that, 
for instance, the matrix C has the elements Cnm)- Lower-case boldface 
letters stand for A^-vectors (so that, for instance, the A^-vector z has the N 
components z„); and the imaginary unit is denoted by i (so that = — 1 , and i 
is of course not a A^-vector!). We occasionally use the Kronecker symbol, with 
its standard definition: dmn = 1 for m = n, 6mn = 0 for m 7 ^ n. And let us 
mention the standard convention according to which an empty sum vanishes 
and an empty product equals unity, i. e. J2f=j ifj) = Oj ^f=J (A) = 1 if 

K <j. m 

A prototypical “goldfish” many-body model [1] [2] [3] [4] is characterized by 
the translation-invariant equations of motion 



A Hamiltonian producing these equations of motion reads as follows: 

Af r N 

^(C;z) = X i+ exp(C«) n {Zn-Zi)~'^ , (lb) 

n=l 1=1, t^n 
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where of course the N coordinates C„ = (t) are the canonical momenta cor¬ 

responding to the canonical particle coordinates Zn = Zn (t). The solution of 
the corresponding initial-values problem is provided by the N roots Zn = Zn (t) 
of the following, rather neat, algebraic equation in the variable z: 


N 

E 

^= 1 , i^n 


Zj (0) -I- i g; Zj (0) 
z - ze (0) 


1 OJ 


exp (i a; t) — 1 


(Ic) 


(Note that this is actually a polynomial equation of degree N in z, as seen by 
multiplying it by nEi - Zrn (0)]). Hence this model is isochronous (when¬ 
ever the parameter oj is positive, as we generally assume hereafter; the special 
case w = 0 is “the” prototypical, nonisochronous, case...): all its solutions are 
completely periodic, with the period T = 27r/a; or, possibly, due to an exchange 
of the particle positions through the motion, with a period that is a (generally 
small, see [5]) integer multiple of T. 

Several solvable generalizations of the goldfish model, characterized by New¬ 
tonian equations of motion featuring additional forces besides those appearing in 
the right-hand side of (la), are known: see for instance [2] [3] [4] and references 
therein. 

Remark 1.1. Above and hereafter we call a many-body model solvable if 
its initial-values problem can be solved by algebraic operations, such as finding 
the N zeros of a known t-dependent polynomial of degree N (of course such an 
algebraic equation can be explicitly solved only for < 4). ■ 

Recently a simple technique has been introduced [ 6 ], which allows to identify 
and investigate additional solvable models of goldfish type; and a few examples 
of such models yielded by this new approach have been identified and tersely 
discussed [ 6 ]. The model treated in this paper is another such new model, which 
is perhaps itself interesting (as all solvable models tend to be), and moreover 
allows—as reported in a separate paper [7]—to obtain remarkable Diophantine 
results for the zeros of (monic) polynomials of degree N the coefficients of which 
are the zeros of Hermite polynomials (see, for instance, [ 8 ]) of degree N. 


2 The model and its solutions 

The Newtonian equations of motion of the new many-body problem of goldfish 
type read as follows: 





N 


n {Zn 
.^=1, 



N 


•E 


/ \N—m 
{Zn) 




N 

i—1^ £^m 


-1 



( 2 a) 
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with 

r m 

c„ = (-ir E n(^.j • ( 2 b) 

1<S1 <S 2 ■ ■ ■ 

Here and hereafter the symbol X)i<si<s 2 <s„<Af signifies the sum from 1 to fV 
over the m (integer) indices sj with j = 1 , ...,m and the restriction si < S 2 < 
... < Sjn', of course this sum vanishes if m > TV, consistently with Notation 
1 . 1 . 

The solutions (t) of this fV-body problem are provided—consistently 

with the expressions (2b)—by the N zeros of the following t-dependent (monic) 
polynomial ippf {z;t) of degree N in z: 

N 

+Y,[c^it) , (3) 

m—1 

where the coefficients Cm (t) are themselves the solutions of the system of N 
ODEs 

N 

Cm — ^ Cm 2 ^ ^ (Cm C^) . (4) 

Because this is a well-known solvable model, the time-dependence of these N 
quantities c™ (t) can be obtained by solving an algebraic (in fact polynomial) 
problem, indeed the solution of the initial-value problem of this dynamical sys¬ 
tem, (4), is provided by the following prescription (see for instance section 4.2.2 
in [4] or [9, 3]): the N quantities Cm = Cm (t) are the N eigenvalues of the NxN 
(t-dependent) matrix 

C(t) = c (0) cos(a; t) + C (0) , (5a) 

LO 

with 

C (0) = diag [cm (0)] , (5b) 

C (0) = diag [c„ (0)] + i [M (0), C (0)] , (5c) 

where of course (see ( 2 b)) 

( m 

Cr^{Q) = {-!)'- E nK(0)] , (5d) 

1<S1<S2...<S™<W ^ 

c,„( 0 ) = (-l)”^ E Eh«,( 0 ) n K( 0 )] , ( 5 e) 

l<Si<S 2 ...< 8 ^<N q=l I ’'= 1 - ’' 5^9 J 

and in the right-hand side of (5c) 

[M (0), C (0)] = M (0) C (0) - C (0) M (0) (5f) 
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with the matrix M (0) defined componentwise in terms of the initial data (0) 
as follows: 


Mnm ( 0 ) = [Zn ( 0 ) - ( 0 )] ^ , U ^ TTl , 

Mnn (0) = - ^ Mni ( 0 ) = - ^ [z„ ( 0 ) - zf ( 0 )]-" . (5g) 

1 =1^ l^n i=l, i^n 

Note that these formulas provide an explicit definition of the N time-dependent 
coefficients Cm (t) in terms of the initial data z„ (0), z„ (0) of the A^-body prob¬ 
lem of goldfish type characterized by the Newtonian equations of motion (2), via 
algebraic operations, amounting essentially to the solution of polynomial equa¬ 
tions of degree N; and that the values z„ (t) at time t of the particle coordinates 
Zn are then provided by the N zeros of the polynomial (z; t), explicitly known 
(see (3)) in terms of its N coefficients Cm (t). It is thereby demonstrated that 
the -/V-body problem of goldfish type characterized by the Newtonian equations 
of motion (2) is solvable (see Remark 1.1). 

Remark 2.1. Let us call attention to a (well known) tricky point associ¬ 
ated with the solution—as described above—of the fV-body problem of goldfish 
type characterized by the Newtonian equations of motion (2). The identifica¬ 
tion of the N eigenvalues of a given matrix is only unique up to permutations, 
and likewise the identification of the zeros of a polynomial is only unique up to 
permutations. Therefore the N coordinates z„ = z„ (t) yielded by the solution 
detailed above are only identified up to permutations of their N labels n. The 
(only) way to identify a specifie coordinate—say, the coordinate zi (<) that cor¬ 
responds to the initial data zi ( 0 ), Zi ( 0 )—is by following the (eontinuous) time 
evolution of the coordinate zi (t) from its (assigned) initial value zi ( 0 ) to its 
value Zi (t) at time t. In this manner one arrives at the uniquely defined value of 
the coordinate zi (t) corresponding to the initial value zi ( 0 ), which coincides of 
course with that uniquely yielded by the time evolution of the A^-body problem 
(2). An analogous phenomenology is also relevant to the solution of model (la). 

■ 

Because of the way system (2) is constructed, its equilibria can be obtained 
by finding the zeros of the polynomials whose eoefficients are equilibria of sys¬ 
tem (4). On the other hand, it is known that the zeros of the A^-th degree 
Hermite polynomial are equilibria of system (4) (see for instance [3]). This 
relationship allows to prove Diophantine properties of the polynomials whose 
coefficients are the zeros of Hermite polynomials. These properties, together 
with the study of the behavior of system (2) in the immediate vicinity of its 
equilibria, are discussed in a separate paper [7]. 

The findings reported above imply the possibility to display in eompletely 
explicit form the solution of the A^-body problem of goldfish type characterized 
by the Newtonian equations of motion (2) for N = 2, 3, 4; but we doubt this 
would be very illuminating, and we therefore leave this task to the eager reader. 
We rather like to emphasize that these findings imply that, for arbitrary N — 
and arbitrary positive w—this A^-body problem is isochronous, all its solutions 
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satisfying the periodicity property 


(t + T) = Zn (t) , (6) 

with T = 2ttI ui (or possibly T might be replaced by a, generally small, integer 
multiple of 27r/a;: see [5]). We end this paper by displaying a few examples of 
this phenomenology. 


Example 2.1. For N = 2 and w = 1, taking into account that ci = —zi — Z 2 
and C2 = 21^2, we see that system (2) reduces to 


2ziZ2 

Zi =- 

2:1 - 2:2 
2iii:2 

Z2 =- 

Zl - Z 2 


1 


{zi - Z2) 
1 


- 


{zi - Z2) 


2(^1 - 1 ) 

(zi + Z2 + ZiZ2)^_ 

_ 2 (z 2 - 1 ) 

(zi + Z2+ ZiZ2)^_ 


(7) 


The assignment w = 1 is motivated by the fact that this quantity, w, can be 
eliminated from system (la) by a constant rescaling of the dependent variables 
Zn(t) and of the independent variable t. 

System (7) has 4 equilibrium configurations {zi\z 2 ^), j = 1,2,3,4, up to 
''( 7 ) • ^( 7 ) 

the exchange of z^ with Z 2 \ whose approximate numerical values are given 
below: 

(i^, 41)) = (0.353553 - 0.762959 i, 0.353553 + 0.762959 i), 

= (-0.54455 + 1.00281 i, 0.54455 - 0.295704 i), , . 

= (-0.54455 - 1.00281 1,0.54455 + 0.295704 i), 

{z[^\z!^'^^) = (-1.26575,0.558645). 

These equilibria can be obtained either by the substitution z„(t) = n = 1,2 
into system (7) and the subsequent solution of the resulting system of algebraic 
equations for zi, Z 2 , or by finding the zeros of the polynomials whose coefficients 
are the equilibria of system (4) for = 2. We note that, in this case where 
N = 2 and w = 1, system (4) has the four equilibria (^)~^)) 

(^,—^), and Two of them are the zeros of the Hermite 

polynomial 772 (c) = 4c^ — 2, which is consistent with the known fact that the 
zeros of the Wth order Hermite polynomial are the equilibria of system (4). 

In Figure 1 each equilibrium of system (7) is represented by the two points 
z^'^ and Z 2 \ labeled by (j), where j = 1,2, 3,4. 

Below we provide graphs of the real and imaginary parts of the components 
of the solution {zi (t) ,Z 2 {t)) of system (7) as functions of time, and some tra¬ 
jectories of the particles zi and Z 2 in the complex z-plane. These graphs have 
been obtained by solving system (7) numerically using Mathematica. 
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Figure 1: Equilibria of system (7). Each equilibrium is represented by the two 
points and 2 : 2 '^ labeled by the index j = 1 , 2 ,3,4, see ( 8 ). 

(7a) Initial conditions zi(0) = + 0.01, ^ 2 ( 0 ) = + 0.01, 

ii(0) = 0.01, i 2 ( 0 ) = -0.01. See Figures 2, 3, 4, 5. 

(7b) Initial conditions Zi(0) = z^"’ + 0.01 + 0.01 i, ^ 2 ( 0 ) = z^ 2 ^ + 0.01 + 0.01 i, 
ii( 0 ) = -0.01, i 2 ( 0 ) = 0.01. See Figures 6 , 7. 

(7c) Initial conditions zi(0) = + 0.01 + 0.01 i, Z 2 ( 0 ) = z^'^ + 0.01 + 0.01 i, 

zi(0) = —0.01, Z 2 ( 0 ) = 0.01. See Figures 8 , 9. 

(7d) Initial conditions Zi(0) = z^^ + 0.01 + 0.01 i, Z2(0) = z^^^ + 0.01 + 0.01 i, 
zi(0) = —0.005, Z 2 ( 0 ) = 0.005. See Figures 10 and 11. 

(7e) Initial conditions Zi(0) = z^^^ + 0.01 + 0.01 i, Z2(0) = z^^^ + 0.01 + 0.01 i, 
zi(0) = —0.01, Z 2 ( 0 ) = 0.01. See Figures 10 and 11. 

Example 2.2. For TV = 3 and w = 1, taking into account that ci = 

—Zi — Z 2 — Z 3 , C 2 = Z 1 Z 2 + Z 1 Z 3 +Z 2 Z 3 , and C 3 = —Z 1 Z 2 Z 3 , we see that system (2) 
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reduces to 


2ziZ2 , 2 ziZ3 

^1 —- 1 - 

Zi - Z2 Zi- 23 

- (^ - z){z - Z ) +-21-F2(-21,2:2,23) +-F3(-21,Z2,23))] , 

2zxZ2 2z2Zz 

Z2 — - 1 - 

Zl - Z2 2^2 - 2:3 

+ 7 - 77 - 7 [zlFi{zi, Z2, z^) + Z2F2 {zi, Z2, z^) + Fz{zi, Z2, z^i))] , 

[Zi - Z2){Z2 - Zz) ^ 

2zxZ-i 2 X 2^3 

Z 3 =- 

■21 - .23 22 - 23 

~7-77-7 [23.^1(21, Z2, 23) + Z3i^2(2l, 22, Z3) +-^3(^1, Z2, 23))] , 

(21 - 23)(Z2 - 23) 


where 


fl(zi, Z2, 23) = Zi + Z2 + 23 - 


(21 + 22 + Z3 — ZiZ223)^ 


(Zi + Z2 + Z3 + Z1Z2 + Z1Z3 + Z223)3 ’ 

^2(21, Z2, 23) = -Z1Z2 - Z1Z3 - Z2Z3 

2 2 
(Zi + Z2 + Z3 + Z1Z2 + Z1Z3 + Z223)3 (ziZ 2 + 21Z3 + Z2Z3 + ZiZ223)3 ’ 
^3(21,22, 23) = Z1Z2Z3 

+ ... 3 )“ ■ (3.33 + .,.3+1^, + J.=3.3)’- 


We obtained equilibria of system (9) as follows. First, we found the equilib¬ 


ria of system (4) for N = 3; they are given by (0, —•y' |) and (0, iy —iy 5 ), 

up to the permutations of the three coordinates. Second, we found the zeros 


of the monic polynomials whose coefficients are the equilibria of system (4) for 
N = S. These zeros are equilibriumum solutions of system (2) because of how 
this system is constructed. Therefore, system (9) has at least 12 equilibrium 
conhgurations (zp\ z^'^^ Zg'^^), j = 1,2,, 12 , up to the permutations of z[^\ 
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''(i) ''(?) 

Z 2 and £3 , whose approximate numerical values are given below: 


{z[^\ = (0.720239 - 0.575751 1,0.720239 + 0.575751 1,-1.44048), 

{zi'^\zi^\zi^^) = (0.397225 + 1.07661 1,-1.12106- 0.854451 1,0.72384- 0.222154 1), 
{zf\z‘i\zf'>) = (0.397225 - 1.07661 1,0.72384 + 0.222154 1,-1.12106 + 0.854451 1), 
{zt\z^2\zf^) = (0.709,-0.3545 - 1.2656 1, -0.3545 + 1.2656 1), 

{zf\z‘i'\zf'>) = (-0.781352,1.00305 - 0.749241 1,1.00305 + 0.749241 1), 
{zf'\zf\zf^) = (0,0.612372 - 0.921816 1,0.612372 + 0.921816 1), 

(z[’^\ zi’'\ zP) = (-0.82853- 0.22063 1,0.82853 - 0.22063 1,1.666 1), 

= (0,-0.673004 + 1.52228 1,0.673004 - 0.297537 1), 

{z[^\zi^\zi^^) = (0.82853 + 0.22063 1,-1.666 1,-0.82853 + 0.22063 1), 

4^°^) = (0, -0.673004 - 1.52228 1, 0.673004 + 0.297537 1), 
(4“\52“\4“^) = (-1.00305 + 0.749241 1,-1.00305 - 0.749241 1,0.781352), 

= (0,-1.87718,0.652438). 

It is possible that the system of algebraic equations characterizing the equilibria 
of (9) has additional solutions besides those listed above. A direct attempt to 
solve this system of algebraic equations using Mathematica was unsuccessful, 
and we did not deem the matter surficiently relevant to justify further investi¬ 
gations. 

In Figure 14 each equilibrium of system (9) is represented by the three points 
z[^\ z^'^ and z^^\ where j = 1,2,..., 12. We leave it to the interested reader to 
compare Figure 14 with the list of equilibria (10), in order to locate the triples 
that correspond to each of the 12 equilibria. 

Below we provide some trajectories^ in the complex plane, of the particles zi, 
Z 2 and Z 3 whose evolution is described by system (9) with the initial conditions 

zi( 0 ) = z®+ 0 . 1 ,ii( 0 )= 0 . 1 , 

Z 2 ( 0 ) = Z® + 0 . 1 ,i 2 ( 0 )= 0 . 1 , 

Z 3 ( 0 ) = 2 f + 0 . 1 ,i 3 ( 0 ) = 0 . 1 , ( 11 ) 

see Figures 15, 16 and 17. These graphs have been obtained by solving sys¬ 
tem (9) with the initial conditions (11) numerically using Mathematica. 




Figure 2: System (7), initial conditions (7a). Graphs of the real and imaginary 
parts of the coordinate zi{t). 



Figure 3: System (7), initial conditions (7a). Trajectory, in the complex z-plane, 
of zi{t). The two dots on the figure indicate the positions of z[^'^ respectively 
2 : 1 ( 0 ), i.e. of a nearby equilibrium point respectively the initial value of Zi{t). 



Figure 4: System (7), initial conditions (7a). Graphs of the real and the imagi¬ 
nary parts of the coordinate 22 (t)- 
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Figure 5: System (7), initial conditions (7a). Trajectory, in the complex z-plane, 
of Z 2 {t). The two dots on the figure indicate the positions of respectively 
Z 2 ( 0 ), i.e of a nearby equilibrium point respectively the initial value of Z 2 {t). 



Figure 6: System (7), initial conditions (7b). Trajectory, in the complex z-plane, 
of Zi{t). The two dots on the figure indicate the positions of z^^^ respectively 
Zi(0), i.e. of a nearby equilibrium point respectively the initial value of Zi{t). 
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Figure 7: System (7), initial conditions (7b). Trajectory, in the complex z-plane, 
of Z 2 {t). The two dots on the figure indicate the positions of respectively 
2 : 2 ( 0 ), i.e. of a nearby equilibrium point respectively the initial value of 22 (t)- 



Figure 8: System (7), initial conditions (7c). Trajectory, in the complex z-plane, 
of Zi{t). The two dots on the figure indicate the positions of z^^ respectively 
21 ( 0 ), i.e. of a nearby equilibrium point respectively the initial value of Zi{t). 
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Figure 9: System (7), initial conditions (7c). Trajectory, in the complex z-plane, 
of Z 2 {t). The two dots on the figure indicate the positions of respectively 
02 ( 0 ), i.e. of a nearby equilibrium point respectively the initial value of 22 (i)- 


3 Outlook 

The interest of the many-body model introduced and discussed above is demon¬ 
strated by its solvable and isochronous character as well as by the remarkable 
trajectories it features - already in the simple N = 2 and N = 3 cases, as shown 
by the graphs reported above. It is also amusing to observe—as the cognoscienti 
will have noted—that this solvable model has been obtained by appropriately 
combining the two solvable equations—see (la) and (4)—which are the pro¬ 
totypes of two, quite different, basic families of solvable many-body problems 
of Newtonian type. Moreover—as reported in [7]—the findings reported above 
provide the point of departure to obtain Diophantine properties of the N ze¬ 
ros of each of the N\ (monic) polynomials the coefficients of which are the N 
zeros of the Hermite polynomial Hn (z) of degree N (these N\ polynomials cor¬ 
respond of course to the N\ permutations of the N zeros of Hn (^)). And a 
rather ample vista of further developments is provided by these findings. 
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Figure 10: System (7), initial conditions (7d). Trajectory, in the complex 2- 

/qX 

plane, of zi{t). The two dots on the figure indicate the positions of z{ and 
2:1(0), i.e. of a nearby equilibrium point respectively the initial value of Zi{t). 



Figure 11: System (7), initial conditions (7d). Trajectory, in the complex 2:- 

/qX 

plane, of 2:2 (t). The two dots on the figure indicate the positions of £2 and 
2:2(0), i.e. of a nearby equilibrium point respecttively the initial value of 2:2 (t). 
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Figure 12: System (7), initial conditions (7e). Trajectory, in the complex z- 
plane, of Zi{t). The two dots on the figure indicate the positions of and 
Zi(0), i.e. of a nearby equilibrium point respectively the initial value of Zi{t). 



Figure 13: System (7), initial conditions (7e). Trajectory, in the complex z- 
plane, of Z2(t). The two dots on the figure indicate the positions of z^"^^ and 
Z2(0), i.e. of a nearby equilibrium point respectively the initial value of Z2(t). 
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Figure 14: Equilibria of system (9). Each equilibrium is represented by the 
three points z[^\ and zjf'^ indexed by j = 1, 2 ,..., 12, see (10). 



Figure 15: System (9), initial conditions (11). Trajectory, in the complex 2- 
plane, of zi (t) . The two dots on the figure indicate the positions of z) ' respec¬ 
tively 21(0), i.e. of a nearby equilibrium point respectively the initial value of 

Zi(t). 



Figure 16: System (9), initial conditions (11). Trajectory, in the complex z- 
plane, of Z 2 {t). The two dots on the figure indicate the positions of and 
Z2(0), i.e. of a nearby equilibrium point respectively the initial value of Z 2 {t). 



Figure 17: System (9), initial conditions (11). Trajectory, in the complex z- 
plane, of Z3(t). The two dots on the figure indicate the positions of z^'^ and 
Z3(0), i.e. of a nearby equilibrium point respectively the initial value of zsit). 
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